Entropy, entanglement, and area: 
analytical results for harmonic lattice systems 
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We revisit the question of the relation between entanglement, entropy, and area for harmonic lattice Hamil- 
tonians corresponding to discrete versions of real free Klein-Gordon fields. For the ground state of the d- 
dimensional cubic harmonic lattice we establish a strict relationship between the surface area of a distinguished 
hypercube and the degree of entanglement between the hypercube and the rest of the lattice analytically, with- 
out resorting to numerical means. We outline extensions of these results to longer ranged interactions, finite 
temperatures and for classical correlations in classical harmonic lattice systems. These findings further suggest 
that the tools of quantum information science may help in establishing results in quantum field theory that were 
previously less accessible. 
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Imagine a distinguished geometrical region of a discretized 
free quantum Klein-Gordon field: what is the entropy asso- 
ciated with a pure state obtained by tracing over the field 
variables outside the region? How does this entropy relate 
to properties of the region, such as volume and boundary 
area? This innocent-looking question is a long-standing is- 
sue indeed, studied in the literature under the key word of 
geometric entropy. Analytical steps supplemented by numeri- 
cal computations for half-spaces and spherical configurations 
in seminal works by Bombelli et al. ll[ and Srednicki |2] 
strongly suggested a direct connection between entropy and 
area. The interest in this quantity for quantum field theory is 
drawn from the fact that geometric entropy is thought to be 
the leading quantum correction to the Bekenstein-Hawking 
black hole entropy |3]. Subsequent work employed vari- 
ous approaches, such as methods from conformal field the- 
ory M, analysis of entropy subadditivity |5] or mode count- 
ing 1 6]. Recently, there has been renewed interest in study- 
ing entanglement and correlations in quantum many-body 
systems and quantum field theory, largely due to availabil- 
ity of novel powerful methods from the quantitative theory 
of entanglement in the context of quantum information the- 
ory |3 11 II 03 [H d Q. Such ideas have previously 
been employed to assess the entanglement in settings of one- 
dimensional spin (see, e.g., Refs. fl2l [l3lo and harmonic 
chains 01 [H. 

This letter gives an analytical answer to the question of 
the scaling of the degree of entanglement for harmonic lat- 
tice Hamiltonians such as discrete versions of the free scalar 
Klein-Gordon field, in arbitrary spatial dimensions. Although 
we encounter a highly correlated system we nevertheless find 
an 'area-dependence' of the degree of entanglement. Our 
analysis is based on methods that have been developed in re- 
cent years in quantum information theory, in particular those 
relating to entanglement in Gaussian (quasi-free) states (see, 
e.g., Ref. 1 14]). These methods allow us to give an analytical 
answer to the question of the scaling of the degree of entan- 
glement between a region and its exterior for harmonic lat- 



tice Hamiltonians such as discrete versions of the free scalar 
Klein-Gordon field, in arbitrary spatial dimensions. It is re- 
markable that although we encounter a highly correlated sys- 
tem, we nevertheless find an "area dependence" of the degree 
of entanglement. 

The Hamiltonian. — The starting point of the argument is a 
discrete lattice version of a free real scalar quantum field. For 
any d > 1 we consider a d-dimensional simple cubic lattice 
n xd comprising n d oscillators. We may write the Hamiltonian 
as 

H = pp T /2 + xVx T /2, (1) 

where x = (xi, . . . , x n d) and p = (px, . . . ,p n d) denote the 
canonical coordinates of the system. The n d x n d -matrix V, 
the potential matrix, specifies the coupling between the oscil- 
lators in the position coordinates. 




FIG. 1 : The harmonic lattice in d — 2 with a distinguished m x m 
region in an n x n lattice. 

For now V will be chosen such that in the continuum 
limit one obtains the Hamiltonian of the real Klein-Gordon 
field, under periodic boundary conditions. We will there- 
fore consider the harmonic lattice Hamiltonian with nearest- 
neighbor interaction. Note that our argument can be extended 
to other types of interactions. The case of next-to-nearest- 



2 



neighbor coupling will also be discussed later in this paper, 
see Ref. IU7I1 for a more general discussion. 

We write V = circ(v) for the circulant matrix whose first 
row is given by the n-tupel v, and also for a block circulant 
matrix where the first block column is specified by a tupel of 
matrices. So in d = 1, we have Vi = circ(l, — c, 0, . . . , 0, — c) 
and in higher dimensions we have a recursive, block-circulant 
structure reflecting rows, layers etc.: 

Vd = circ(Vd-i, — cl n d-i,0, ... ,0, — c\ n d~i) 

with a < 2cd < 1. From now on we will write V instead of 
Vd- 

Entanglement and area dependence. — We denote the 
ground state of the system by p. For a distinguished cubic 
region m xd in a lattice n xd (see Fig.Q its entropy of entan- 
glement is 



E n ,r. 



tr p„, m l0g p n ,r 



The reduced density matrix p„ iTO is formed by tracing out the 
variables outside the region m xd . We will show the following: 

The entropy of entanglement of the distinguished region m xd 
in the lattice n xd satisfies 



lim En 



Q(m d - 1 ), 



(2) 



where is the Landau-theta. More specifically, we have that 
Cim d ~ 1 < E n . m < C2 m' 1 ^ 1 for sufficiently large m, with 
appropriate C\ , C'2 > 0. 

The 'area dependence' manifests itself as follows: For a 
linear chain, the entropy of entanglement is bounded by quan- 
tities that are independent of the size of the distinguished in- 
terval. In two dimensions, this dependence is linear in the 
length of the boundary, in three dimensions to the area of the 
boundary. Indeed, one can show that while all oscillators are 
correlated with all oscillators, the correlations over the bound- 
ary decay very quickly. In effect, for fixed interaction strength 
11911 . the only significant contribution comes from within a 
finite width, the correlation length, along the boundary, and 
thus leads to a surface dependence of the correlations. This 
intuition forms the basis of the following, fully analytical ar- 
gument, where the above statement is proven by finding upper 
and lower bounds for which the statement holds. 

The upper bound. — The ground state p of the coupled 
harmonic system in Eq. {0 is a Gaussian (quasi-free) state 
with vanishing first moments. The second moments of p can 
be collected in the covariance matrix 7, which is defined as 
7 i:fe = 2Re[RjR k p] for j,k = l,...,2n d , where R = 
(xi, . . . , x n d,pi, . . . ,p„d) is the vector of canonical coordi- 
nates. In terms of the potential matrix V the covariance ma- 
trix of the ground state is then found to be 7 = V~ x / 2 © V x l 2 
1121 • From entanglement theory we know that an upper bound 
for the entropy of entanglement is provided by the logarithmic 
negativity Em — ln||p r || tr , where p r is the partial transpose 
of p, and || • || t , denotes the trace norm fil. Following Ref. [9] 
we find 



E 



N 



■max(0,Aj(Q- 1))), 



where Xj (Q) are the non-increasingly ordered eigenvalues of 
the matrix 



y-l/2p V l/2p 



(4) 



In a reordered list of canonical coordinates (such that the inner 
oscillators are counted first) P is the diagonal matrix P = 
— l m d ® t n d_ m d and the potential matrices can be written as 



A B 
B T C 



V 1 ' 2 



D E 
E T F 



T 



E 
E T 



The matrices B and E describe the couplings between the m d 
oscillators forming the distinguished hypercube and the rest 
of the lattice. On using V~ 1 I 2 V 1 I 2 = 1, we arrive at 

Q-l = -2V~ 1/2 T 

This is convenient as it will turn out that the detailed structure 
of V~ x l 2 will not have to be considered and we can concen- 
trate on the properties of the matrix T. To avoid taking the 
maximum in Eq. Q we bound the eigenvalues by their abso- 
lute values, 



+ ^((2-1)1) 



<^|A,-(Q-l)| = ||Q-l|| tr , 

i=i 

where we have employed that ln(l + x) < x for all x > 0. 
Since the trace norm is unitarily invariant 1 16], we may further 
write 



\\Q - l|| tr - 2 \\V-^ 2 T\\ lr < 2 AifV- 1 ") ||T|| tr . 

Here we also have that V^ 1 ^ 2 is symmetric. The spectrum of 
V can be obtained via discrete Fourier transform and yields 
XiiV- 1 ' 2 ) - (1 -2cd)- 1 /2. 

Now the trace norm of T can be bounded from above by 
the sum of the absolute values of all the matrix elements of T, 
which is known as the l\ matrix norm 1 16]. Therefore, 



l/2x 



En < 



= y 

2cd .4< 



T 



In the following we will bound the matrix elements of V 1 ^ 2 
and consequently those of T. The explicit implementation 
of the multidimensional discrete Fourier transform is non- 
technical yet involved. To achieve a more compact nota- 
tion, we introduce the lattice coordinate vectors k, 1 where 
kj, lj = 0, . . . , n— 1 and j = 0, . . . , d— 1. For the considered 
lattice structure we may write Vk 1 = V^d-i . i ^d-i , , 

for the interaction term between site k and 1. The matrix ele- 
ments of V 1 / 2 are then given by 



(3) V. 



1/2 



k.l 



^-1 27Tifc' (fc,- -/,)/« 

k' j=0 



^ 2nk< 
2c > cos 



1/2 



r=0 



n 
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To bound these, we replace the square root by its power 
series expansion in the parameter 2c. This converges if 
2cd < 1, which coincides with the constraint imposed by 
the positivity of the potential matrix. We will use (1 — 
z)V2 = _ <p« i b s x s , with < B 8 < 1 and the fact that 
12q=i e 2 * lpq l n — for integer p and q unless p is a multi- 
ple of n. With this the non-diagonal elements of V 1 / 2 , and 
analogously V^ 1 ^ 2 , are bounded by 



y s(k,l) 



> v. 



-1/2 



k,l 



> 0> K 



1/2 



k.l 



> -'- 



^(k.l) 



(5) 



where s(k, 1) = (ko — lo) H h (fc n -i ^ ^n-i), y = 2cg? and 

< kj — lj < n/2. This demonstrates the exponential decay 
of the off-diagonal elements in these block circulant matrices. 
The remaining matrix elements are determined by the periodic 
boundary conditions under the exchange kj — lj i— * n — (kj — 
lj ) . Note that generally s (k, 1) is simply the number of lattice 
steps one has to make starting at site k to reach site 1. If, for 
example, s(k, 1) = 1 then the oscillators are direct neighbors. 

We may now proceed with the computation of the l\ norm 
of T, i.e., of the blocks in V 1 / 2 that describe the coupling 
between the distinguished region and the rest of the lattice. 
Given that the region is a hypercube, this can be done in a 
transparent way. Consider the set Co of m d — (to — 2) d oscil- 
lators of the hypercube that lie directly on the boundary and 
successively the sets C r of (to — 2r) d — (to — 2r — 2) d oscil- 
lators inside that are exactly r steps away from the surface of 
the hypercube. Starting from the set Co and taking s steps on 
the lattice one can reach less than (m + 2s) d — m d oscillators 
outside the hypercube m xd . Therefore we find that the sum 
of all the elements of T that couple oscillators from the set Co 
to oscillators outside the hypercube is bounded by 



So < 2^2{{m + 2s) d — TO 



l 



y 



Now consider the contribution from the set C%. Clearly, any 
oscillator outside the hypercube that can be reached from Cf- 
in s + k steps can be reached from £ m s steps. Therefore 
we can bound the sum Sk of all the elements of T that couple 
the set Ck to oscillators outside the hypercube by 



S k < 2 J2 ((m + 2(s-k)) d - 



s=k+l 



As a consequence we obtain 



m/2 



E N < 



< 



k=Q 



m/2 



£((TO + 2 S ) d _ m rf)JL_£y 



Using the binomial expansion of (m + 2s) d and the Gamma- 
function to bound expressions of the form X)^o 2/ s (2s) fe we 
find for to > Ad/\ ln(y)| the bound 



En < 



IQd 



VI - 2cd(l - 2cd) 2 \ ln(l - 2cd)\- 



m d -\ (6) 



which is the desired upper bound that is linear in the number 
of oscillators on the surface of the hypercube. 

Lower bound. — In the following we demonstrate that the 
degree of entanglement, measured by the entropy of entan- 
glement, is asymptotically at least linear in the number of 
oscillators. The entropy of entanglement depends only on 
the symplectic spectrum of the covariance matrix 7^4 corre- 
sponding to the reduced Gaussian state of the interior. The 
non-increasingly ordered symplectic eigenvalues satisfy = 
(Ai(AD)) 1 / 2 > 1 from which the entropy of entanglement 
can be evaluated as 



Pi + l . Ht + l fii - 1 . Mi - 1 



log 



log; 



For Hi > 1 each bracketed terms in the sum can be bounded 
from below by log /ij. Because /Ltj < ((1 + 2c) /(l — 2c)) 1 / 4 
for all i, we find 

m d , ra d 

S > £ log (1 + ( Mi - 1)) > £ - 1) . 

i— 1 — 1 



Employing that for = (yl + k — l)/fc we have a/T-FIc > 
1 + /3x in x € [0, k] and A, (AD) = 1 + \{-BE T ) we find 



X^-BE^) /ii-1 V 7 

The factors in front of the trace can be bounded from above 
by a quantity that is independent of both to and n. All the ele- 
ments of V~ 1/>2 and of — E are positive. Using the techniques 
that led to Eqs. 0) we find 

, ±V2| > i /cyO" 1 * _JL_ 

1 k -' 1 " 2 V2/ 1 -c 2 



As a consequence we have 



k.l 



Now we employ counting methods analogous to those used 
in the derivation for the upper bound we find an expression 
linear in the area. We take into account only contributions to 
the above sum that correspond to the 2dm d ~ 1 oscillators that 
can be reached in each step moving outwards orthogonal to 
the surface of the hypercube. We thus obtain a lower bound 
proportional to the surface of the hypercube m xd for to > too 
and appropriate m . This concludes the proof. 

In the following we will briefly describe possible extensions 
of the above results that can be obtained by similar techniques, 
including more general interactions, thermal states and classi- 
cal correlations in classical systems. 

'Squared interactions'. — The basic intuition behind the 
entanglement-area dependence becomes most transparent for 
the specific class of interactions for which the potential matri- 
ces V is of the form V — W 2 with a circulant band-matrix 
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W. In that case the covariance matrix of the ground state is 
given by 7 = W -1 © W. In this case one arrives at Eq. 
since one can show that (i) the number of terms contributing 
to the symplectic spectrum of the reduced covariance matrix 
is linear in the number of degrees of freedom at the bound- 
ary of the region, and (ii) the respective symplectic eigen- 
values are bounded from above and below independently of 
n and m. Note that property (i) is equivalent to the exis- 
tence of a 'disentangling' symplectic unitary transformation 
local to inside and outside of the regions such that only os- 
cillators near to the boundary remain entangled. Taking e.g. 
Vi = circ(l + 2c 2 , -2c, c 2 , 0, . . . , 0, c 2 , -2c) - the case of 
nearest-neighbor and smaller next-to-nearest-neighbor inter- 
actions - allows to show that only the oscillators exactly at 
the boundary contribute to the logarithmic negativity and that 
Ai(Q) < 2/(1 - 2c) - 1, with Q being defined as in Eq. 0. 
For the same interaction in d = 1 spatial dimension one can 
even exactly calculate the symplectic spectrum of the reduced 
covariance matrix by means of a simple recursion relation. 
In the limit m — ► 00 this results in the two non-vanishing 
symplectic eigenvalues fix = fi 2 = (1 — c 2 /^ 2 ) -1 / 2 , where 
q = c+ 1/2 ± (c+ 1/4) 1 / 2 . 

Entanglement and area in classical systems. — It should 
be noted that, perhaps surprisingly, an 'area-dependence' 
can also be established analytically for classical correlations 
in classical harmonic lattice systems 11711. It is notewor- 
thy that this result on classical systems can be established 
most economically using quantum techniques namely, map- 
ping the problem onto that of a quantum harmonic lattice with 
a squared interaction as has been described above. 

Entanglement and area at finite temperature. — The prop- 
erty of squared interactions leading to effective disentangle- 
ment extends to thermal states and therefore permits the proof 
of the linear entanglement-area dependence for finite temper- 



atures. In that case operational entanglement measures such 
as the distillable entanglement have to be used. They can 
be bounded from below by the hashing inequality and above 
again by the logarithmic negativity 11711 . 

Summary and outlook. — For certain harmonic lattice 
Hamiltonians, e.g. discrete versions of the real Klein-Gordon 
field, we have proven analytically that the degree of entan- 
glement between a hypercube and its environment can be 
bounded from above and below by expressions proportional to 
the number of degrees of freedom on the surface of the hyper- 
cube . This establishes rigorously a connection between en- 
tanglement and area in this system. Intuitively, this originates 
from the fact that one can approximately decouple the oscilla- 
tors in the interior and the exterior up to a band of the width of 
the order of the correlation length of the system, which can be, 
as outlined for the case of next-to-nearest neighbor coupling, 
equal to just one lattice unit. 

Our results can be extended to a wide variety of harmonic 
lattice Hamiltonians, both quantum and classical, and a future 
publication |17] will present details for more general inter- 
actions, both ground and thermal states and a careful discus- 
sion of the continuum limit, where the effective interaction 
strength is modified. These results in particular rely in an 
essential way on the insights and techniques that have been 
obtained in recent years in the development of a quantitative 
theory of entanglement in quantum information science. 
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